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ABSTRACT ; 

* 

A simple dynamic finite element algorithm for analyzing a propagating 
mixed mode crack tip Is presented. A double nodlng technique# which can be 
easily incorporated Into existing dynamic finite element codes# Is used 

A 

together with a corrected J Integral to extract modes I and II dynamic stress 
Intensity factors of a propagating crack. Tha utility of the procedure Is 
demonstrated by analyzing test problems Involving a mode I central crack 
propagating In a plate subjected to uniaxial tension# a mixed mode I and II 
stationary# slanted central crack In a plate subjected to uniaxial Impact 
loading# and a mixed mode I and II extending# slanted single edge crack In a 
plate subjected to uniaxial tension. 
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INTRODUCTION 

Most of the recent numerical studies In dynamic fracture mechanics are 
restricted to mode I crack propagation with the crack extending along the 
line of symmetry. Crack extension under such condition can be modeled with 
finite element method by simply releasing the crack tip node at the line of 
symmetry following a prescribed nodal force versus crack extension history. 
Mode I dynamic stress Intensity factor In dynamic fracture analysis of 
Isotropic material Is Justified since the crack will ordinarily propagate In 
the direction perpendicular to the maximum principal stress direction. 

There are conditions* however* when the crack will deviate from its 
self-similar crack extension path to curve Cl] or to bifurcate [2] under 
specific static or dynamic loadings. Such crack extension away from the line 
of symmetry In a finite element mesh cannot be accomplished by the above- 
mentioned simple nodal release mechanism. The double-nodlng C3] and the 
nodal -grafting [4] techniques are two procedures which have been used to 
model off-axis crack propagation. While details of the former are not 
available* the latter requires a higher order isoparametric element In the 
dynamic finite element code. The double nodlng technique presented In this 
paper was developed for use with an Implicit dynamic finite element code 
which utilizes a conventional Isoparametric quadrilateral element. 

THEORETICAL BACKGROUND 
Double Nodlng Technique 

_ Consider a slanted crack In a two dimensional space with a local x-y 
coordinate system oriented along the local crack tip region as shown In 
Figure 1. The dynamic equations of motion for the crack tip element with 
nodal displacements of {q}* nodal velocities of (q)* nodal accelerations of 
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{q}» nodal forces of (F), stiffness of CK3 and mass matrix of CM3 are: 

CK3(q) + CM] (q) » {FJ (1) 

If the 1th and jth degrees of freedom in the above displacements* velocltl 
and accelerations are constrained to be equal 

q i = v <1 * v q i * q j <2> 

through double nodlng* then define th6 relative displacement {q * > ?s 

{q} = [T] {q"}, etc. (3) 

{F") = [T] T {F) <4) 

where T = 1 when m = n 
mn 

= 1 when m = 1 and m » j 
= 0 all other m and n 


and 


q "l “ q 1 ~ q J “ 0 When k * 1 


* % 


when k f- 1 


( 5 ) 


Substituting Equations <3>* (4) and (5) Into Equation (1) yields 
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» {F} 

= [T] T (F) 

■ [T] T [K]{q) + [T] T [M]«) 

= [T] T tK][T]{q'} + [T] T [M][T]{q'> (6) 


The equivalent stiffness and nass matrices In the above equivalent dynamic 
equations of motion are 


[K'] = [T] T [K][T] 

(7) 

[in - [T] T [M][T] 

(8) 


Equations (3) through (8) constitutes the equivalent system for the slanted 
crack with two and more degrees of freedom constrained by the double noded 
current and future crack tip nodes. The equivalent nodal displacements are 
determined by applying Newmark's beta method to the dynamic equations of mo- 
tion. The nodal displacements at the double nodes are then determined by 
Equation (3). 

Nodal Force Release 

Crack extension with a traction free crack surface requires elimination 
of the nodal forces normal# F 2 and F 2p , and tangential# Fj^ and F lp » to the 
crack surface as shown In Figure 1. The crack tip nodal force components 
are released simultaneously in equal increments to model a linear translation 
of the crack tip to its new crack tip node# (Fj p # Fj^» F 2p » F^} ■ tF lr’ V 
F 2p # F 2 £>(1- ^ Kwhere bar denotes the forces at the double node prior to 
crack extension# A is the current crack tip location and d Is the crack 
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distance to Its adjacent node. Details of this release procedure as well as 
Its comparison with others are described In Reference C53. 


. Integral 

A 

The corrected J Integral was used to determine the modes I and II 

A 

dynamic stress intensity factors# Kj and Kjj# respectively. The J Integral 
as defined by Kishimoto et al C63 Is 



p Vu dA 


(9) 


where the Indices refer to the local and X2 coordinates shown In Figure 1# 
W Is the strain energy density# n^ Is the surface normal to the Integration 
path and P Is the density. 


By taking a symmetric integration path with respect to the crack 
extension direction# as shown In Figure 1# the J Integral# which contains 

A 

both Kj and Kjj# can be decomposed Into two path-independent Integrals Jj and 
J "n •* 

A A A 

J = Jj + J n (10) 


A A * 

where Jj and Jjj are only functions of Kj and respectively. Detailed 
formulations of Jj and Jjj are described by Klshlmoto et al. £63. 

For a stationary crack# Ishlkawa et al [73 has shown that 


£l K 2 
Bi r K i 


(11a) 
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^ K+l " v 2- 

J II ‘ 8p ~U~ 


• 

(lib) 

where . 

3 - 4v 

i 

plane strain 


k = : i--‘ 

3-v 

T+v 

plane stress 



and V and v are the shear modulus and Poisson's ratio# respectively. 
For a propagating crack# Atlurl et al. C8»9] has shown that 


. Fj(a) 2 
Jr = “4 Kf 

I 2 \i I 


(12a) 


" 2 
l = y c 

'll 2y *11 


(12b) 


where Fj (a) 


SiO-So) , % 0+S,} 2 

- L jjr L - 


2*2 


-2(S,-S 2 ). ( ,+S 2* 


/(HS^O+Sz) 


(12c) 


F„(a) ■ 


s 2 (l-s 2 ) (1+S.) . 2 

-hr 1 - t 2S 2< 1+s i>-2sr-° +s 2> 


(HS 2 )" 


-2(S 2 -S 1 )-==4= 

2 1 /(1+S,)(1+S 2 ) (12d) 
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S 


2 

1 



2 


1 ‘ (^) 2 


(12a) 


D = 4S^2 - 


(1+S 2 ) 2 


(12f) 


and a , C._ and are the crack velocity, dllatational and dlstortlonal stress 
wave velocities, respectively. 

NUMERICAL PROCEDURE 

The numerical procedure consists of inputting the above double nodlng 
technique to an ii? -licit dynamic finite element code and releasing the crack 
tip nodal force 11nt>-rly In accordance with a prescribed crack tip motion. 

As shown by Equations (12), the associated dynamic stress Intensity factors, 

A A 

Kj and K^, can he determined by numerically evaluating the Jj and J II 
Integrals along predetermined symmetric contours surrounding the Instantane- 

A 

ous crack tip. Separation of the J Integral of Equation (10) Into the two 

A A 

J , and Integrals Is accomplished by using the decomposition procedure 
Initially developed by Ishikawa et al. C73 for a static mixed mode crack and 
which was later extended to a dynamically loaded mixed mode crack by Klshlmo- 
to et al. C6D. The procedure consists of first decomposing the displace- 
ments, strains, stresses and body forces Into symmetric and anti -symmetric 

A A 

components after which Equation (9) Is used to compute Jj and Jjj Integrals, 
respectively. In the following, three problems which utilize the above 

A A 

mentioned double nodlng technique together with the Jj and 3^ integral 
technique of determining Kj and Kjj are described. The examples Involved 
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steel plate under a plane strain state of stress tilth material properties of 
V - 2.94 x 10 10 N/m 2 , v - 0.286 and P - 2.45 x 10 3 Kg/e 3 . 

EXTENDING CENTRAL CRACK IN A PLATE SUBJECTED TO UNIAXIAL TENSION 

The Broberg problem [103 of a central crack extending at constant 
velocity In a uniaxial tension field has been used by many to verify their 
dynamic finite element codes. This problem yields only a mode I dynamic 
stress Intensity factor* but It offers the only comparison between results 
generated by other dynamic finite element codes. Figure 2 shows the finite 
element break down of a half-plane as well as the seven Integration paths 
used in this analysis. The complete specimen geometry Is shown In the 
legend of Figure 3. Variations In the numerically determined mode I dynamic 
stress Intensity factors with the Integration paths for five crack tip 
locations of a/W * 0.2, 0.3# 0.4, 0.5 and C.6 and three crack velocities of 
a/C 2 = 0, 0.2 and 0.6 are shown In Table 1, where a Is the half crack length 
and W Is the half plate width. Maximum difference of 1 percent between the 
Kj computed for the various Integration paths shows that the path 
Independency Is satisfied for all practltal purpose. 

Figure 3 shows the changes In the static and the dynamic stress 'nten- 
slty factors for two dynamic crack velocities. Also shown for comparison 
purpose are the static results of Islda C113 and the two dynamic results of 
Nlshloka and Atlurl Cl2,133. The 13 percent lower Kj at the shorter crack 
length of a/W » 0.25 at a higher crack velocity of a/C 2 ■ 0.6 Is due In part 
to the coarser and conventional element used In this analysis. Otherwise, 
good agreements between the various results are notd. 


STATIONARY SLANTED CENTRAL CRACK IN A PLATE SUBJECTED TO UNIAXIAL IMPACT 
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LOADING 

Mixed mode dynamic stress Intensity factors In a stationary slanted 
central crack In a plate subjected to uniaxial Impact loading of oH(t) was 
first analyzed by Thau and Lu C14] using Wlener-Hopf technique and by 
Klshlmoto et al. [6] using finite element analysis. This problem was also 
analyzed in this paper In order to compare the applicability of present 
finite element algorithm with the above known results. Figure 4 shows the 
finite element break down and the five Integration paths used In this 
analysis. Table 2 shows the variations in the normalized K^ and Kjj at four 
time intervals of the stationary crack. While the differences In K n for the 
five Integration paths Is as high as 8 percent# the relatively small Kjj 
values makes this difference insignificant. Figure 5 shows changes in tran- 
sient Kj and Kjj with time for the stationary cr»rk where R Q » R^# R 2 denoting 
the first arrival times at the crack tip for Incident wave from leading edge# 
the reflected waves from upper/lower boundaries and the reflected wave from 
the rear edge# respectively. While agreements between the three results for 
the stationary crack Is good at the Initial stage, the finite element results 
differ towards the latter part of time Interval of 14 ^ 20 s. Fqually 
puzzling Is the large differences between the static stress Intensity factors 
for the plate under static loading by these two finite element analyses. 

Such discrepancy may result from different finite element mesn and geometry 
(triangular versus quadrilateral) by Klshlmoto et al. C6] and by the present 
computation. Nevertheless# the general results of these two analyses are 
very similar. 


EXTENDING SLANTED SINGLE EDGE CRACK IN A PLATE SUBJECTED TO UNIAXIAL TENSION 
As a further study In mixed mode dynamic crack# an extending slanted 
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single edge crack In a plate subjected to uniaxial tension* which Is the 
dynamic counterpart of the static solution of Bowie C153 » was analyzed. This 
crack was extended along Its original crack direction at two crack velocities 
of a/C-. * 0.2 and 0.6. Figure 6 shows the finite element breakdown as well 
as the four Integration paths used In this analysis. Table 3 shows the 
variations In Kj and Kjj at five crack tip locations for the static and two 
propagating cracks. Unlike the previous case, little variation In Kjj are 
noted, possibly due to the relatively larger values of Kjj In this problem. 
Figure 7 shows the changes In Kj and Kjj with crack extension. Also shown 
for comparison Is the static results of Bowie [153 which are within 4X of the 
present analysis. While one would not expect self-similar cracx propagation 
under this loading condition. It is Interesting to note the closeness of Kj 
and Kjj values at the higher crack velocity of a/C 2 = 0.6. 

CONCLUSIONS 

A double nodlng technique suitable for dynamic finite element analysis 
of a crack extending under mixed mode loading condition Is presented. Mixed 

A 

mode dynamic stress Intensity factors are computed by using the J Integrals 
as modified by Atlurl et al. [8,93 for the extending crack. 

The procedure was used to determine the dynamic stress intensity factors 
of an extending central crack In a unlaxlilly loaded plate, the dynamic 
stress Intensity factors of a stationary slanted central crack In an unlaxl- 
ally Impacted plate and the static and dynamic stress Intensity factors of an 
extending slanted, single edge crack In a uniaxial ly loaded plate. The 
computed values by this procedure were generally In good agreement with known 
results. 

When modified by appropriate criteria of dynamic crack curving C13 and 
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branching C2]# this simple procedure can be used to determine the dynamic 
fracture parameters associated with such problems. 
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TABLE 3^ 

PATH INDEPENDENCE OF 3 INTEGRAL FOR 
EXTENDING SLANTED SINGLE EDGE CRACK PROBLEM 


Kj/a/TT* Kjj/cf/Sw 




Path 




Path 




a/C 2 

a/W 

1 

2 

3 

4 

1 

2 

3 

4 


C.2 

0.4314 

0.4437 



0.2243 

0.2220 




0.3 

0.6583 

0.6638 

0.6726 

0.6691 

0.3226 

0.3212 

0.3198 

0.3270 

0 

0.4 

1.0079 

1 7111 

1.0167 

1.0245 

0.4498 

0.4485 

0.4455 

0.4405 


0.5 

1.5784 

1.5829 

1.5897 

1.5972 

0.6134 

0.6102 

0,6036 

0.5865 


0.6 

2.5244 

2.5451 

2.5226 

2.5223 

0.8462 

0.8230 

0.7749 

0.7488 


0.3 

0.4205 

0.4252 

0.4316 

0.4254 

0.2708 

0.2708 

0.2708 

0.2749 


0.4 

0.5497 

0 5550 

0.5563 

0.5613 

0.3122 

0.3111 

0.3099 

0.3108 

0.2 

0.5 

0.7549 

0. ‘5 90 

0.7583 

0.7627 

0.4163 

0.4168 

0.4167 

0.4142 


0.6 

0.8831 

O.^JlO 

0.8902 

0.9001 

0.5579 

0.5387 

0.5227 

0.5263 


0.3 

0.1762 

0.1863 

0.1914 

0.1774 

0.2166 

0.2178 

0.2161 

0.2252 

0.6 

0.4 

0.2095 

0.2240 

0.2187 

0.2190 

0.2287 

0.2322 

0.2360 

0.2326 


0.5 

0.2211 

0.2347 

0.2321 

0.2351 

0.2661 

0.2648 

0.2641 

0.2639 


0.6 

0.2571 

0.2793 

0.2750 

0.2759 

0.2282 

0.2834 

0.2772 

0.2783 


a is the crack length. 

W is the plate width. 
a is the applied uniaxial load 
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FIGURE I. DOUBLE NODING, NODAL FORCE RELEASE * 
AND J INTEGRAL PATH, T. 








FIGURE 2. FINITE ELEMENT BREAKDOWN OF A HALF 
PLANE OF AN EXTENDING CENTRAL CRACK. 
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FIGURE 3. NORMALIZED STRESS INTENSITY FACTORS OF 
EXTENDING CENTRAL CRACK. 
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FIGURE 4. FINITE ELEMENT BREAKDOWN OF STATIONARY 
SLANTED CENTRAL CRACK. 
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SLANTED CENTRAL CRACK. 
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FINITE ELEMENT BREAKDOWN OF EXTENDING SLANTED 
SINGLE EDGE CRACK. 
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FIGURE 7. NORMALIZED STRESS INTENSITY FACTORS OF 
EXTENDING SLANTED SINGLE EDGE CRACK. 



